Comprehensive analyses of a large human gut Bacteroidales culture collection reveal species and strain level diversity and evolution.

Species of the Bacteroidales order are among the most abundant and stable bacterial members of the human gut microbiome with diverse impacts on human health. While Bacteroidales strains and species are genomically and functionally diverse, order-wide comparative analyses are lacking. We cultured and sequenced the genomes of 408 Bacteroidales isolates from healthy human donors representing nine genera and 35 species and performed comparative genomic, gene-specific, mobile gene, and metabolomic analyses. Families, genera, and species could be grouped based on many distinctive features. However, we also show extensive DNA transfer between diverse families, allowing for shared traits and strain evolution. Inter- and intra-specific diversity is also apparent in the metabolomic profiling studies. This highly characterized and diverse Bacteroidales culture collection with strain-resolved genomic and metabolomic analyses can serve as a resource to facilitate informed selection of strains for microbiome reconstitution.

Dietary- and host-derived metabolites are used by diverse gut bacteria for anaerobic respiration.

Respiratory reductases enable microorganisms to use molecules present in anaerobic ecosystems as energy-generating respiratory electron acceptors. Here we identify three taxonomically distinct families of human gut bacteria (Burkholderiaceae, Eggerthellaceae and Erysipelotrichaceae) that encode large arsenals of tens to hundreds of respiratory-like reductases per genome. Screening species from each family (Sutterella wadsworthensis, Eggerthella lenta and Holdemania filiformis), we discover 22 metabolites used as respiratory electron acceptors in a species-specific manner. Identified reactions transform multiple classes of dietary- and host-derived metabolites, including bioactive molecules resveratrol and itaconate. Products of identified respiratory metabolisms highlight poorly characterized compounds, such as the itaconate-derived 2-methylsuccinate. Reductase substrate profiling defines enzyme-substrate pairs and reveals a complex picture of reductase evolution, providing evidence that reductases with specificities for related cinnamate substrates independently emerged at least four times. These studies thus establish an exceptionally versatile form of anaerobic respiration that directly links microbial energy metabolism to the gut metabolome.

Activity of Gut-Derived Nisin-like Lantibiotics against Human Gut Pathogens and Commensals.

Recent advances in sequencing techniques unveiled the vast potential of ribosomally synthesized and post-translationally modified peptides (RiPPs) encoded in microbiomes. Class I lantibiotics such as nisin A, widely used as a food preservative, have been investigated for their efficacy in killing pathogens. However, the impact of nisin and nisin-like class I lantibiotics on commensal bacteria residing in the human gut remains unclear. Here, we report six gut-derived class I lantibiotics that are close homologues of nisin, four of which are novel. We applied an improved lantibiotic expression platform to produce and purify these lantibiotics for antimicrobial assays. We determined their minimal inhibitory concentration (MIC) against both Gram-positive human pathogens and gut commensals and profiled the lantibiotic resistance genes in these pathogens and commensals. Structure-activity relationship (SAR) studies with analogs revealed key regions and residues that impact their antimicrobial properties. Our characterization and SAR studies of nisin-like lantibiotics against both pathogens and human gut commensals could shed light on the future development of lantibiotic-based therapeutics and food preservatives.

Fecal metabolite profiling identifies liver transplant recipients at risk for postoperative infection.

Metabolites produced by the intestinal microbiome modulate mucosal immune defenses and optimize epithelial barrier function. Intestinal dysbiosis, including loss of intestinal microbiome diversity and expansion of antibiotic-resistant pathobionts, is accompanied by changes in fecal metabolite concentrations and increased incidence of systemic infection. Laboratory tests that quantify intestinal dysbiosis, however, have yet to be incorporated into clinical practice. We quantified fecal metabolites in 107 patients undergoing liver transplantation (LT) and correlated these with fecal microbiome compositions, pathobiont expansion, and postoperative infections. Consistent with experimental studies implicating microbiome-derived metabolites with host-mediated antimicrobial defenses, reduced fecal concentrations of short- and branched-chain fatty acids, secondary bile acids, and tryptophan metabolites correlate with compositional microbiome dysbiosis in LT patients and the relative risk of postoperative infection. Our findings demonstrate that fecal metabolite profiling can identify LT patients at increased risk of postoperative infection and may provide guideposts for microbiome-targeted therapies.

An evolution-based framework for describing human gut bacteria.

The human gut microbiome contains many bacterial strains of the same species ('strain-level variants'). Describing strains in a biologically meaningful way rather than purely taxonomically is an important goal but challenging due to the genetic complexity of strain-level variation. Here, we measured patterns of co-evolution across >7,000 strains spanning the bacterial tree-of-life. Using these patterns as a prior for studying hundreds of gut commensal strains that we isolated, sequenced, and metabolically profiled revealed widespread structure beneath the phylogenetic level of species. Defining strains by their co-evolutionary signatures enabled predicting their metabolic phenotypes and engineering consortia from strain genome content alone. Our findings demonstrate a biologically relevant organization to strain-level variation and motivate a new schema for describing bacterial strains based on their evolutionary history.

Bifidobacteria metabolize lactulose to optimize gut metabolites and prevent systemic infection in patients with liver disease.

Progression of chronic liver disease is precipitated by hepatocyte loss, inflammation and fibrosis. This process results in the loss of critical hepatic functions, increasing morbidity and the risk of infection. Medical interventions that treat complications of hepatic failure, including antibiotic administration for systemic infections and lactulose treatment for hepatic encephalopathy, can impact gut microbiome composition and metabolite production. Here, using shotgun metagenomic sequencing and targeted metabolomic analyses on 847 faecal samples from 262 patients with acute or chronic liver disease, we demonstrate that patients hospitalized for liver disease have reduced microbiome diversity and a paucity of bioactive metabolites, including short-chain fatty acids and bile acid derivatives, that impact immune defences and epithelial barrier integrity. We find that patients treated with the orally administered but non-absorbable disaccharide lactulose have increased densities of intestinal bifidobacteria and reduced incidence of systemic infections and mortality. Bifidobacteria metabolize lactulose, produce high concentrations of acetate and acidify the gut lumen in humans and mice, which, in combination, can reduce the growth of antibiotic-resistant bacteria such as vancomycin-resistant Enterococcus faecium in vitro. Our studies suggest that lactulose and bifidobacteria serve as a synbiotic to reduce rates of infection in patients with severe liver disease.

A ubiquitous mobile genetic element disarms a bacterial antagonist of the gut microbiota.

DNA transfer is ubiquitous in the gut microbiota, especially among species of Bacteroidales. In silico analyses have revealed hundreds of mobile genetic elements shared between these species, yet little is known about the phenotypes they encode, their effects on fitness, or pleiotropic consequences for the recipient's genome. Here, we show that acquisition of a ubiquitous integrative and conjugative element encoding an antagonistic system shuts down the native contact-dependent antagonistic system of Bacteroides fragilis . Despite inactivating the native antagonism system, mobile element acquisition increases fitness of the B. fragilis transconjugant over its progenitor by arming it with a new weapon. This DNA transfer causes the strain to change allegiances so that it no longer targets ecosystem members containing the same element yet is armed for communal defense.

Virulence and genomic diversity among clinical isolates of ST1 (BI/NAP1/027) Clostridioides difficile.

Clostridioides difficile produces toxins that damage the colonic epithelium, causing colitis. Variation in disease severity is poorly understood and has been attributed to host factors and virulence differences between C. difficile strains. We test 23 epidemic ST1 C. difficile clinical isolates for their virulence in mice. All isolates encode a complete Tcd pathogenicity locus and achieve similar colonization densities. However, disease severity varies from lethal to avirulent infections. Genomic analysis of avirulent isolates reveals a 69-bp deletion in the cdtR gene, which encodes a response regulator for binary toxin expression. Deleting the 69-bp sequence in virulent R20291 strain renders it avirulent in mice with reduced toxin gene transcription. Our study demonstrates that a natural deletion within cdtR attenuates virulence in the epidemic ST1 C. difficile isolates without reducing colonization and persistence. Distinguishing strains on the basis of cdtR may enhance the specificity of diagnostic tests for C. difficile colitis.

Virulence and genomic diversity among clinical isolates of ST1 (BI/NAP1/027) Clostridioides difficile.

Clostridioides difficile (C. difficile) , a leading cause of nosocomial infection, produces toxins that damage the colonic epithelium and results in colitis that varies from mild to fulminant. Variation in disease severity is poorly understood and has been attributed to host factors (age, immune competence and intestinal microbiome composition) and/or virulence differences between C. difficile strains, with some, such as the epidemic BI/NAP1/027 (MLST1) strain, being associated with greater virulence. We tested 23 MLST1(ST1) C. difficile clinical isolates for virulence in antibiotic-treated C57BL/6 mice. All isolates encoded a complete Tcd pathogenicity locus and achieved similar colonization densities in mice. Disease severity varied, however, with 5 isolates causing lethal infections, 16 isolates causing a range of moderate infections and 2 isolates resulting in no detectable disease. The avirulent ST1 isolates did not cause disease in highly susceptible Myd88 -/- or germ-free mice. Genomic analysis of the avirulent isolates revealed a 69 base-pair deletion in the N-terminus of the cdtR gene, which encodes a response regulator for binary toxin (CDT) expression. Genetic deletion of the 69 base-pair cdtR sequence in the highly virulent ST1 R20291 C. difficile strain rendered it avirulent and reduced toxin gene transcription in cecal contents. Our study demonstrates that a natural deletion within cdtR attenuates virulence in the epidemic ST1 C. difficile strain without reducing colonization and persistence in the gut. Distinguishing strains on the basis of cdtR may enhance the specificity of diagnostic tests for C. difficile colitis.

Immunomodulatory fecal metabolites are associated with mortality in COVID-19 patients with respiratory failure.

Respiratory failure and mortality from COVID-19 result from virus- and inflammation-induced lung tissue damage. The intestinal microbiome and associated metabolites are implicated in immune responses to respiratory viral infections, however their impact on progression of severe COVID-19 remains unclear. We prospectively enrolled 71 patients with COVID-19 associated critical illness, collected fecal specimens within 3 days of medical intensive care unit admission, defined microbiome compositions by shotgun metagenomic sequencing, and quantified microbiota-derived metabolites (NCT #04552834). Of the 71 patients, 39 survived and 32 died. Mortality was associated with increased representation of Proteobacteria in the fecal microbiota and decreased concentrations of fecal secondary bile acids and desaminotyrosine (DAT). A microbiome metabolic profile (MMP) that accounts for fecal secondary bile acids and desaminotyrosine concentrations was independently associated with progression of respiratory failure leading to mechanical ventilation. Our findings demonstrate that fecal microbiota composition and microbiota-derived metabolite concentrations can predict the trajectory of respiratory function and death in patients with severe SARS-Cov-2 infection and suggest that the gut-lung axis plays an important role in the recovery from COVID-19.

KHSRP loss increases neuronal growth and synaptic transmission and alters memory consolidation through RNA stabilization.

The KH-type splicing regulatory protein (KHSRP) is an RNA-binding protein linked to decay of mRNAs with AU-rich elements. KHSRP was previously shown to destabilize Gap43 mRNA and decrease neurite growth in cultured embryonic neurons. Here, we have tested functions of KHSRP in vivo. We find upregulation of 1460 mRNAs in neocortex of adult Khsrp-/- mice, of which 527 bind to KHSRP with high specificity. These KHSRP targets are involved in pathways for neuronal morphology, axon guidance, neurotransmission and long-term memory. Khsrp-/- mice show increased axon growth and dendritic spine density in vivo. Neuronal cultures from Khsrp-/- mice show increased axon and dendrite growth and elevated KHSRP-target mRNAs, including subcellularly localized mRNAs. Furthermore, neuron-specific knockout of Khsrp confirms these are from neuron-intrinsic roles of KHSRP. Consistent with this, neurons in the hippocampus and infralimbic cortex of Khsrp-/- mice show elevations in frequency of miniature excitatory postsynaptic currents. The Khsrp-/- mice have deficits in trace conditioning and attention set-shifting tasks compared Khsrp+/+ mice, indicating impaired prefrontal- and hippocampal-dependent memory consolidation with loss of KHSRP. Overall, these results indicate that deletion of KHSRP impairs neuronal development resulting in alterations in neuronal morphology and function by changing post-transcriptional control of neuronal gene expression.

Transcriptional changes and preservation of bone mass in hibernating black bears.

Physical inactivity leads to losses of bone mass and strength in most mammalian species. In contrast, hibernating bears show no bone loss over the prolonged periods (4-6 months) of immobility during winter, which suggests that they have adaptive mechanisms to preserve bone mass. To identify transcriptional changes that underlie molecular mechanisms preventing disuse osteoporosis, we conducted a large-scale gene expression screening in the trabecular bone and bone marrow, comparing hibernating and summer active bears through sequencing of the transcriptome. Gene set enrichment analysis showed a coordinated down-regulation of genes involved in bone resorption, osteoclast differentiation and signaling, and apoptosis during hibernation. These findings are consistent with previous histological findings and likely contribute to the preservation of bone during the immobility of hibernation. In contrast, no significant enrichment indicating directional changes in gene expression was detected in the gene sets of bone formation and osteoblast signaling in hibernating bears. Additionally, we revealed significant and coordinated transcriptional induction of gene sets involved in aerobic energy production including fatty acid beta oxidation, tricarboxylic acid cycle, oxidative phosphorylation, and mitochondrial metabolism. Mitochondrial oxidation was likely up-regulated by transcriptionally induced AMPK/PGC1α pathway, an upstream stimulator of mitochondrial function.

Three Strains of Tobacco etch virus Distinctly Alter the Transcriptome of Apical Stem Tissue in Capsicum annuum during Infection.

Tobacco etch virus (TEV; genus Potyvirus) is flexuous rod shaped with a single molecule of single-stranded RNA and causes serious yield losses in species in the Solanaceae. Three TEV strains (HAT, Mex21, and N) are genetically distinct and cause different disease symptoms in plants. Here, a transcriptomic RNA sequencing approach was taken for each TEV strain to evaluate gene expression of the apical stem segment of pepper plants during two stages of disease development. Distinct profiles of Differentially Expressed Genes (DEGs) were identified for each TEV strain. DEG numbers increased with degree of symptom severity: 24 from HAT, 1190 from Mex21, and 4010 from N. At 7 days post-inoculation (dpi), when systemic symptoms were similar, there were few DEGs for HAT- and Mex21-infected plants, whereas N-infected plants had 2516 DEGs. DEG patterns from 7 to 14 dpi corresponded to severity of disease symptoms: milder disease with smaller DEG changes for HAT and Mex21 and severe disease with larger DEG changes for N. Strikingly, in each of these comparisons, there are very few overlapping DEGs among the TEV strains, including no overlapping DEGs between all three strains at 7 or 14 dpi.

Mfd Affects Global Transcription and the Physiology of Stressed Bacillus subtilis Cells.

For several decades, Mfd has been studied as the bacterial transcription-coupled repair factor. However, recent observations indicate that this factor influences cell functions beyond DNA repair. Our lab recently described a role for Mfd in disulfide stress that was independent of its function in nucleotide excision repair and base excision repair. Because reports showed that Mfd influenced transcription of single genes, we investigated the global differences in transcription in wild-type and mfd mutant growth-limited cells in the presence and absence of diamide. Surprisingly, we found 1,997 genes differentially expressed in Mfd- cells in the absence of diamide. Using gene knockouts, we investigated the effect of genetic interactions between Mfd and the genes in its regulon on the response to disulfide stress. Interestingly, we found that Mfd interactions were complex and identified additive, epistatic, and suppressor effects in the response to disulfide stress. Pathway enrichment analysis of our RNASeq assay indicated that major biological functions, including translation, endospore formation, pyrimidine metabolism, and motility, were affected by the loss of Mfd. Further, our RNASeq findings correlated with phenotypic changes in growth in minimal media, motility, and sensitivity to antibiotics that target the cell envelope, transcription, and DNA replication. Our results suggest that Mfd has profound effects on the modulation of the transcriptome and on bacterial physiology, particularly in cells experiencing nutritional and oxidative stress.

Identification of regenerative processes in neonatal spinal cord injury in the opossum (Monodelphis domestica): A transcriptomic study.

This study investigates the response to spinal cord injury in the gray short-tailed opossum (Monodelphis domestica). In opossums spinal injury early in development results in spontaneous axon growth through the injury, but this regenerative potential diminishes with maturity until it is lost entirely. The mechanisms underlying this regeneration remain unknown. RNA sequencing was used to identify differential gene expression in regenerating (SCI at postnatal Day 7, P7SCI) and nonregenerating (SCI at Day 28, P28SCI) cords +1d, +3d, and +7d after complete spinal transection, compared to age-matched controls. Genes showing significant differential expression (log2FC ≥ 1, Padj ≤ 0.05) were used for downstream analysis. Across all time-points 233 genes altered expression after P7SCI, and 472 genes altered expression after P28SCI. One hundred and forty-seven genes altered expression in both injury ages (63% of P7SCI data set). The majority of changes were gene upregulations. Gene ontology overrepresentation analysis in P7SCI gene-sets showed significant overrepresentations only in immune-associated categories, while P28SCI gene-sets showed overrepresentations in these same immune categories, along with other categories such as "cell proliferation," "cell adhesion," and "apoptosis." Cell-type-association analysis suggested that, regardless of injury age, injury-associated gene transcripts were most strongly associated with microglia and endothelial cells, with strikingly fewer astrocyte, oligodendrocyte and neuron-related genes, the notable exception being a cluster of mostly downregulated oligodendrocyte-associated genes in the P7SCI + 7d gene-set. Our findings demonstrate a more complex transcriptomic response in nonregenerating cords, suggesting a strong influence of non-neuronal cells in the outcome after injury and providing the largest survey yet of the transcriptomic changes occurring after SCI in this model.

Corrigendum: Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis.

[This corrects the article DOI: 10.3389/fpls.2016.00762.].

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5' untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5' UTR. As expected, mutations in double-stranded (DS) regions of the 5' UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5' cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Experimental drought reduces genetic diversity in the grassland foundation species Bouteloua eriopoda.

Understanding the resistance and resilience of foundation plant species to climate change is a critical issue because the loss of these species would fundamentally reshape communities and ecosystem processes. High levels of population genetic diversity may buffer foundation species against climate disruptions, but the strong selective pressures associated with climatic shifts may also rapidly reduce such diversity. We characterized genetic diversity and its responsiveness to experimental drought in the foundation plant, black grama grass (Bouteloua eriopoda), which dominates many western North American grasslands and shrublands. Previous studies suggested that in arid ecosystems, black grama reproduces largely asexually via stolons, and thus is likely to have low genetic variability, which might limit its potential to respond to climate disruptions. Using genotyping-by-sequencing, we demonstrated unexpectedly high genetic variability among black grama plants in a 1 ha site within the Sevilleta National Wildlife Refuge in central New Mexico, suggesting some level of sexual reproduction. Three years of experimental, growing season drought reduced black grama survival and biomass (the latter by 96%), with clear genetic differentiation (higher FST) between plants succumbing to drought and those remaining alive. Reduced genetic variability in the surviving plants in drought plots indicated that the experimental drought had forced black grama populations through selection bottlenecks. These results suggest that foundation grass species, such as black grama, may experience rapid evolutionary change if future climates include more severe droughts.

Cranberry-derived proanthocyanidins induce a differential transcriptomic response within Candida albicans urinary biofilms.

Candida albicans is one of the most common causes of hospital-acquired urinary tract infections (UTIs). However, azoles are poorly active against biofilms, echinocandins do not achieve clinically useful urinary concentrations, and amphotericin B exhibits severe toxicities. Thus, novel strategies are needed to prevent Candida UTIs, which are often associated with urinary catheter biofilms. We previously demonstrated that cranberry-derived proanthocyanidins (PACs) prevent C. albicans biofilm formation in an in vitro urinary model. To elucidate functional pathways unique to urinary biofilm development and PAC inhibition, we investigated the transcriptome of C. albicans in artificial urine (AU), with and without PACs. C. albicans biofilm and planktonic cells were cultivated with or without PACs. Genome-wide expression analysis was performed by RNA sequencing. Differentially expressed genes were determined using DESeq2 software; pathway analysis was performed using Cytoscape. Approximately 2,341 of 6,444 total genes were significantly expressed in biofilm relative to planktonic cells. Functional pathway analysis revealed that genes involved in filamentation, adhesion, drug response and transport were up-regulated in urinary biofilms. Genes involved in carbon and nitrogen metabolism and nutrient response were down-regulated. In PAC-treated urinary biofilms compared to untreated control biofilms, 557 of 6,444 genes had significant changes in gene expression. Genes downregulated in PAC-treated biofilms were implicated in iron starvation and adhesion pathways. Although urinary biofilms share key features with biofilms formed in other environments, many genes are uniquely expressed in urinary biofilms. Cranberry-derived PACs interfere with the expression of iron acquisition and adhesion genes within urinary biofilms.

Histone Citrullination Represses MicroRNA Expression, Resulting in Increased Oncogene mRNAs in Somatolactotrope Cells.

Peptidylarginine deiminase (PAD) enzymes convert histone arginine residues into citrulline to modulate chromatin organization and gene expression. Although PADs are expressed in anterior pituitary gland cells, their functional role and expression in pituitary adenomas are unknown. To begin to address these issues, we first examined normal human pituitaries and pituitary adenomas and found that PAD2, PAD4, and citrullinated histones are highest in prolactinomas and somatoprolactinomas. In the somatoprolactinoma-derived GH3 cell line, PADs citrullinate histone H3, which is attenuated by a pan-PAD inhibitor. RNA sequencing and chromatin immunoprecipitation (ChIP) studies show that the expression of microRNAs (miRNAs) let-7c-2, 23b, and 29c is suppressed by histone citrullination. Our studies demonstrate that these miRNAs directly target the mRNA of the oncogenes encoding HMGA, insulin-like growth factor 1 (IGF-1), and N-MYC, which are highly implicated in human prolactinoma/somatoprolactinoma pathogenesis. Our results are the first to define a direct role for PAD-catalyzed histone citrullination in miRNA expression, which may underlie the etiology of prolactinoma and somatoprolactinoma tumors through regulation of oncogene expression.